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ABSTRACT 


Presented in this paper was an overview on research works on solar radiation basics and photovoltaic 
generation. Also, a complete PV modeling and investigation on the effect of using multi-axes sun-tracking 
systems on the electrical generation was carried out to evaluate its performance using the case study of 
the Monastir city, Tunisia. The effects of azimuth and tilt angles on the output power of a photovoltaic 
module were investigated. The instantaneous increments of the output power generated by a photo- 
voltaic module mounted on a single and dual-axis tracking system relative to a traditional fixed panel 
were estimated. The results show that the yearly optimal tilt angle of a fixed panel faced due to the south 
is close to 0.9 times Monastir latitude. The gain made by the module mounted on a single-axis tracking 
panel relative to a traditional fixed panel was analyzed. The monthly increments of the gain are more 
noticeable for two critical periods which correspond to those surrounding the summer and the winter 


Gain solstice dates. It reaches the value of 10.34% and 15% in the summer and winter solstice periods, respec- 
tively. However, the yearly gain relative to a fixed panel installed with the yearly optimal tilt angle is 
5.76%. In some applications, covering loads at early morning or late afternoon hours and in order to more 
optimize the solar systems exploitation suggest the adjustment of the PV panel orientation to azimuth 
angles different from the south direction by using a dual-axis tracking installation. The gain made by 
this recommendation relative to a traditional fixed panel is evaluated. This gain reaches 30% and 44% 

respectively in the winter and summer solstice days. 
© 2011 Elsevier Ltd. All rights reserved. 
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1. Introduction 


Renewable energy resources have enormous potential and can 
satisfy the present world energy demand. They can enhance diver- 
sity in energy supply markets, secure long-term sustainable energy 
supply, and reduce local and global atmospheric emissions [1]. One 
of the important renewable energy is the solar energy which can 
be captured anywhere and can be directly converted into electric 
power through PV panels [2] which allow the production of elec- 
tricity using photovoltaic cells that convert solar irradiation into 
useable electricity and, as a result, have become one of the most 
reliable sources of low pollutant energy. Besides, PV energy produc- 
tion represents an environmentally advantageous and sustainable 
method of maintaining an energy intensive standard of living [3-5]. 

The use of photovoltaic panels was originally planned as a type 
of technology to be used for external space applications [6-11]. 
The PV applications have been harnessed in various areas such as 
water pumping, rural electrification, remote meteorological sta- 
tions, maritime and railway crossing, production of hydrogen, grid 
connection systems, refrigeration, low power rural industrial appli- 
cations, etc. [12-26]. 

On the other hand, recent technological advancements have 
stimulated the development and the use of PV panels in the com- 
mercial market for its use in building appliances. Regardless of 
its requisite high primary costs, PV systems have been adopted 
in many countries as an efficient means of producing energy in 
rural and urban areas that show a higher output in rural areas [27]. 
Thereby, the worldwide annual development rate in PV industry 
has averaged more than 30% in the last decade [28,29]. 

APV module is composed of several solar cells. The output power 
of a PV module is mainly based on two factors, i.e. cell temperature 
and solar radiation incident on it [30]. The intensity of incident solar 
radiation on a panel is affected by the installation azimuth and tilt 
angle, as both angles influence the incident angle of sunlight on 
it. The optimal angles vary with conditions such as the geographic 
latitude, climate, atmospheric composition, utilization period and 
so forth [31,32]. Itis known that the optimum value of yearly max- 
imum output energy can be obtained from PV modules oriented 
facing south and the optimum tilt angle is equal to the local lat- 
itude [33,34]. However the optimal tilt angle varies significantly 
with the geographic latitude, climate condition, utilization period 
of time, etc. [35,36]. 

It is widely acknowledged that Southern orientation for PV pan- 
els in the Northern Hemisphere is optimal and because that result 
indicate that an average of 98.6% a system’s performance with the 


optimal angle can be obtained using the latitude angle for the tilted 
panel [37]. 

Many empirical correlations are available in literature for 
estimating the optimum tilt angle under various operational con- 
ditions. It is found that the optimum value of yearly maximum 
output energy can be obtained from PV modules oriented facing 
south with a tilt angle in the range of 20-30° [38]. Duffie and 
Beckman [33] suggested the use of (latitude + 15°) +15° where the 
plus and minus signs correspond to the use in winter and sum- 
mer, respectively. Chow and Chan [39] analyzed the irradiation 
data measured in a coastal region of South China and found that, 
over a year, solar collectors should be placed with a tilt angle of 
2.8° greater than the latitude with the azimuth facing south-west. 
Shariah et al. [40] showed that the yearly optimum inclination 
angle for solar collectors in Jordan is less than the latitude by about 
5-8°. 

Gunerhan and Hepbasli [41] found that the optimal angles for 
solar collectors in Izmir, Turkey were in good agreement with those 
calculated by the expressions proposed by Nijegorodov et al. [42] 
which consist of a 12 expressions set for determining the monthly 
optimal tilt angle for latitudes between 60° north and 60° south. 
However, Yakup and Malik [43] found that the optimum tilt angles 
in Brunei Darussalam were almost less than the ones calculated 
by [42], due to the change in radiation patterns from location to 
location. 

The performance analysis mentioned above are not yet suf- 
ficient, because they use usually monthly mean values of daily 
solar irradiation on inclined surfaces. However, in many cases, 
these values are not sufficient, because they do not allow a pre- 
cise idea of the different energy phenomena, which take place 
in the heart of the production system (inertia phenomenon, 
shadowing masks, etc.). So, it is very important to study the instan- 
taneous performance of PV arrays at different tilt angles and 
azimuths. 

In this paper, an evaluation study and an overview on research 
works of all over the solar basics and photovoltaic generation 
are presented. The mentioned models can be considered as the 
representative of worldwide average. Its details can be found 
therein. Furthermore, many factors that affect the instantaneous 
PV production are analyzed for specific days in different situ- 
ations using the case study of Monastir city. Also, the amount 
of extra output power generated by a single-axis tracked panel 
relative to fixed panel, installed with the yearly optimal tilt 
angle which is computed before, are estimated in this paper. 
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Zenith 


Fig. 1. Zenith angle, slope, surface azimuth angle and solar azimuth angle for a tilted 
surface. 


Finally, the gains made by a dual-axis tracked panel are evalu- 
ated. 


2. Solar radiation basic: an overview 


The sun, with a diameter of 1.39 x 109 m, emit a total energy out- 
put of 3.8 x 102° MW, however only a tiny fraction of 1.7 x 10!4 kW 
is captured by the earth after traveling a mean distance of 
1.495 x 10!! m [44]. The eccentricity of the earth’s orbit is such the 
distance between the sun and the earth varies by 1.7% [33]. Fur- 
ther sun-earth geometric relationships will be illustrated in the 
following section. Also, the direction from which beam radiation is 
received and its angle of incidence on various surface over any time 
span will be analyzed. Then, the last part of the following section is 
provided to describe the modeling used to estimate the total solar 
radiation on horizontal and tilt surface. 


2.1. Sun-earth geometric relationships and incidence angle 


The geometric relationships between a plane of any particular 
orientation relative to the earth at any time and the beam solar 
radiation can be described by appropriate relationships in term of 
several angles [33] illustrated in Fig. 1. 


2.1.1. True and mean anomaly of the earth 

The true anomaly is an angular parameter that defines the posi- 
tion of the earth moving along a Keplerian orbit. It is the angle 
between the direction of the periapsis and the current position of 
the earth as seen from the sun [45]. 

The mean anomaly B, called also the fractional year, is a parame- 
ter relating position and time for the earth moving in a Kepler orbit. 
It is what the true anomaly would be if the earth moved with con- 
stant speed along a perfectly circular orbit around the sun in the 
same time. 

The mean anomaly B is given by the following formula: 


_ 360(n— 1) 


B= — sas (1) 


where n is the ordinal date. n=1 at January 1st and n=365 at 
December 31st. 


2.1.2. Earth obliquity and sun declination 
The earth obliquity, also called earth axial tilt, is the angle 
between its rotational axis, and a perpendicular to the orbital plane. 


The tilt and direction of the rotational axis with respect to the orbit 
change very slowly (over thousands of years) [46]. The actual value 
of obliquity is 23.45° [47]. 

The declination of the sun 6 is the angle between the rays of the 
sun and the plane of the earth’s equator. Solar declination varies 
with the seasons and its period is one year. At the solstices, ô reaches 
its maximum value of 23.45°’. Therefore 6=+23.45° at the sum- 
mer solstice while 6 =—23.45° at the winter summer solstice. At 
the equinoxes, ô is equal to 0°. 

In term of the mean anomaly B, the declination is given by [48]: 


ô = 0.006918 — 0.3999912 cos B + 0.070257 sin B 
— 0.006758 cos 2B + 0.0000907 sin 2B 
— 0.002679 cos 3B + 0.00148 sin 3B (2) 


In term of ordinal date, the declination is expressed by the fol- 
lowing formula [49]: 


ne 284+n 
ô = 23.45° sin (sco2*) (3) 


2.1.3. Equation of time 

The equation of time is the difference between mean time (time 
based on the average motion ofa fictional sun which moves at a uni- 
form rate) and true time (time based on a non-uniform motion of 
the physical sun around the earth). The true sun’s apparent motion 
varies due to the elliptical nature of the earth’s orbit and the incli- 
nation of the axis of the earth’s rotation. It is calculated using the 
following expression: 


E = 3.8210-°(75 + 1868 cos B — 32077 sin B 
— 14619 cos 2B — 40890 sin 2B) (4) 


2.1.4. Latitude and longitude 

The latitude, g, is the angular location north or south of the 
equator, north positive; —90° < g < 90+. 

The longitude À is a geographic coordinate that specifies the 
east-west position of a point on the Earth’s surface. g and A are 
angular measurements, usually expressed in degrees, minutes and 
seconds. 


2.1.5. Solar time and hour angle 

Standard time zone is defined by geometrically subdividing the 
earth’s spheroid into 24lunes, bordered by meridians each 15°. 
Time zones are based on Greenwich Mean Time (GMT), the mean 
solar time at longitude O° (the Prime meridian). Thus, to convert 
civil time, tc, to solar time ts, two corrections should be applied. The 
first one is the correction for the difference in longitude between 
the observer’s meridian and the meridian on which the civil time is 
based. The second correction is from the equation of time (4). Solar 
time is expressed as follows: 


ts te+ ie —Zc+E (5) 
where Zc is the time zone east of GMT. 

The hour angle w is the solar angular displacement east or west 
of the local meridian due to the rotation of the earth (15° per hour). 
The hour angle is negative in the morning and positive afternoon. 
Therefore, the hour angle can be written as follows: 


w = (ts — 12)- 15° (6) 


4056 T. Maatallah et al. / Renewable and Sustainable Energy Reviews 15 (2011) 4053-4066 


2.1.6. Zenith, solar altitude and azimuth angles 

The zenith angle 0; is defined as the angle between the vertical 
and the line to the sun, which is, the angle of incidence of a beam 
radiation on a horizontal plane. The zenith angle is expressed in 
terms of latitude ¢, declination ô and hour angle, w as follows: 


cos 6; = cos ọ cos ô cos w + sin g sin ô (7) 


The solar altitude angle as is the angle between the horizontal 
and the line to the sun, that is, the complement of the zenith angle. 

The solar azimuth angle ys is defined as the angular dis- 
placement from south of the projection of beam radiation on 
horizontal plane, with zero due south, east negative, west positive; 
—180° < y < 180°. 


2.1.7. Surface azimuth and slope angles 

The surface azimuth angle y is the deviation of the projection on 
a horizontal plane of the normal to the surface from the local merid- 
ian. The slope, £, referred also to the surface inclination, is the angle 
between the plane of the surface in question and the horizontal; 
0° < 6 < 180. 


2.1.8. Angle of incidence 

The incidence angle 0 is the angle between the beam radia- 
tion and the normal to the surface with any orientation. It can be 
expressed using the following equation [33]: 


cos 8 = sin g sin ô cos B 
— cos g sin ô sin f cos y 
+ cos Q cos ô COS f cos w (8) 
+ sin ọ cos ô sin f cos y cos w 
+ cos ô sin f sin y sin œw 


2.2. Extraterrestrial solar radiation 


2.2.1. Solar constant 

The solar constant, called also the extraterrestrial radiation, 
Gsc is defined as the energy from the sun per unit time received 
on a unit area of surface perpendicular to the direction of prop- 
agation of the radiation at mean earth-sun distance outside the 
atmosphere. In earlier works, the estimations and measurements 
of the extraterrestrial irradiance have re-examined many times. In 
fact, Johnson revised the Abbot’s extraterrestrial irradiance value of 
1322 W/m? to 1395 W/m? (1954) [50]. In 1971, the NASA and the 
American Society of Testing and Materials accepted the Thekaekara 
and Drummond [51], value of 1353 W/m?. In 1978, Frohlich [52] 
recommended a new value of 1373 W/m2. Wilson et al. reported 
1368 W/m? as a new value in 1982. In this study, the 1367 W/m? 
value adopted by the World Radiation Center (WCR) is used. 


2.2.2. Normal extraterrestrial radiation 

The normal extraterrestrial radiation Gon is defined as the 
amount of solar radiation striking a surface normal to the sun’s 
rays at the top of the earth’s atmosphere. Spencer [48] provided 
an accurate equation of the normal extraterrestrial radiation as 
follows: 


Gon = Gsc(1.00011 + 0.034221 cos B + 1.0001280 sin B 
+0.000719 cos 2B + 0.000077 sin 2B (9) 


The normal extraterrestrial radiation can be given by a simple 
adequate equation for most engineering calculations: 


Gon = Gse (1 + 0.033 cos (10) 


360n 
365 ) 


2.2.3. The horizontal extraterrestrial radiation 

The horizontal extraterrestrial radiation G,, is the amount of 
solar radiation received by a horizontal surface at the top of atmo- 
sphere. The horizontal extraterrestrial radiation can be calculated 
basing on Eqs. (7) and (10) as follows: 


Goh = Gon Cos 0z (11) 


The average extraterrestrial horizontal radiation over a time 
step can be calculated as: 


(@2-@1) . a; 
80 sin osin ô 


(12) 


Goh = = [cos g cos &(sinw, — sinw,)+7 


where œw and w are respectively the hour angle at the beginning 
of the time step and the hour angle at the end of the time step. 


2.3. Radiation on horizontal surface 


2.3.1. Global radiation _ 
The global radiation, G on horizontal plane, is the sum of diffuse 
radiation Gg and the direct radiation Gp. It is expressed as follows: 


G = Gp + Ga (13) 


2.3.2. Diffuse radiation 

The diffuse radiation G, is solar reaching the earth’s surface after 
having been scattered from the direct solar beam by molecules or 
suspensoids in the atmosphere [53]. 

The clearness index factor K is defined as the ratio of a particular 
global radiation to the horizontal extraterrestrial radiation: 


G 
Goh 
Basing on the data from four U.S. and one Australian station 
and using the clearness index K, Erbs et al. [54] (1982) suggested 


a correlation to resolve global radiation into its beam and diffuse 
components: 


1.0 — 0.09K 

(for K < 0.22) 
_ J 0.9511 — 0.1604K + 4.388K? — 16.638K? + 12.336K4 (15) 
— ) (for0.22 < K < 0.8) 

0.165 


(for K > 0.8) 


K = (14) 


OO 


Furthermore, Reindl et al. [55] (1979) proposed another corre- 
lation as follows: 


1.0 — 0.258K 

(for K < 0.3) 

1.45 — 1.67K 

(for 0.3 < K < 0.78) 
0.147 

(for K > 0.78) 


a] 
I 
T 
o 


Moreover, Mondol et al. [56] (2008) suggested the following 
correlation: 


0.98K 

(forK < 0.2) 

0.61092 + 3.6259K — 10.171K2 + 6.388K? 
(for 0.22 < K < 0.8) 

0.672 — 0.474K 

(for K > 0.8) 


OO 
I 
pa 
y 


2.3.3. Beam radiation 

The beam radiation, G, is defined as solar radiation that travels 
from the sun to the earth’s surface without any scattering by the 
atmosphere. In this study, the correlation proposed by Erbs et al. 
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[54] will be used to estimate the diffuse radiation G4 and then the 
beam radiation G; is calculated by subtracting the diffuse radiation 
from the global horizontal radiation. 


2.4. Radiation on inclined surface 


Many authors presented a several models in order to convert 
the solar radiations for a horizontal surface to that for a tilted one. 
These models use the same method of calculating beam and ground 
reflected radiation on a tilted surface and the only difference exists 
in the treatment of the diffuse radiation [53]. 


2.4.1. Beam radiation 
The beam radiation on tilted surface, G,;, can be estimated by 
multiplying its value on a horizontal surface, Gp, by a geometric 
factor Rp which depends on the zenith and incidence angles. Rp is 
given by the following the expression [33,57,58]: 
cos 0 


R, = cos 0, (18) 


Thus, the beam radiation on tilted surface Gp; is expressed as: 


cos 0 


Gpt = Gp x Rp = G 
be = Go x Ro =G X oe, 


(19) 


2.4.2. Diffuse radiation 

As mentioned previously, the only difference among the models 
appears in the assessment of sky-diffuse component [59-64]. 

It has been pointed out that the sky diffuse component is con- 
sidered to be the largest potential source of computational error 
[65,66]. In fact, the most complex problem of estimating the diffuse 
irradiance on tilted surface, Gy, is related to the fact that this dif- 
fuse component does not have an isotropic distribution in the sky 
and this distribution is not uniform over the time. Many authors 
investigated the diffuse radiation component basing on various 
measurements from several stations. In order to evaluate this com- 
ponent, they used many models classified basically on three types: 
isotropic, circum-solar and anisotropic. 


2.4.2.1. Isotropic models. The isotropic models presume that the 
intensity of diffuse radiation is uniform over the entire hemisphere 
sky dome which could be true in the case of overcast skies. Hence, 
the diffuse radiation incident on a tilted surface depends on the 
fraction of the sky dome seen by it [53]. 

Later, Dorota [67] revealed that many models calculating 
solar radiation on an inclined surface of any particular orienta- 
tion have been performed and the isotropic diffuse sky model, 
Hottel-Woertz-Liu-Jordan [68,69], is the most used. According to 
this model diffuse component is given as follows: 


Gat = Ga x (=) (20) 


where £ is the slope of the considered surface. 

However, despite the popularity of the isotropic models, the- 
oretical as well as experimental results have shown that the 
simplified calculations of isotropic model do not give a real picture 
of the isolation conditions [67]. In fact, they showed the poorest 
performance and should not be used for estimating the diffuse 
radiation on a tilted surface [55]. 


2.4.2.2. Circum-solar models. The circum-solar model presumes 
thatall the radiation comes from the direction of the sun and its sur- 
roundings, and, therefore, that the diffuse component is calculated 
analogously to the direct component [33,57]: 


Gat = Ga x Rp (21) 


Skydome Circumsolar diffuse 


Isotropic diffuse 


Horizon brightening diffuse 


Ground reflected 


Fig. 2. Schematic view of the solar radiation distribution on a tilted surface accord- 
ing to anisotropic diffuse sky model. 


Unfortunately, this model could only be applied in the case of 
completely clear skies, generally overestimates the diffuse compo- 
nent of the radiation [57]. 


2.4.2.3. Anisotropic models. The anisotropic model assume the 
anisotropy of the diffuse sky radiation in the circumsolar region 
(sky near the solar disc) and the isotropy diffuse component dis- 
tribution from the rest of the sky dome. This model is widely 
accepted, it is used both for clear and cloudy or partly cloudy days 
[57]. Based on the above, an important number of models consid- 
ering anisotropic distribution of the diffuse irradiance in the sky 
are proposed [70-74]. In this context, one of the most confirmed 
anisotropic sky model is the HDKR, Hay—Davies—Klucher-Reindl, 
model [67]. This model assumes that there are three components 
to the diffuse solar radiation: an isotropic component which comes 
from all parts of the sky equally, a circumsolar diffuse component 
which is concentrated in the sky near the sun and a horizon bright- 
ening component which emanates from the sky near the horizon 
(Fig. 2). 

According to The HDKR model the diffuse component radiation 
incident on a tilted surface is as follows: 


Gat = Ĝa(RA:)+ (1 a (ge) Í +f sin? (5) (22) 


where the factor f accounts the fact that more diffuse radiation 
comes from the horizon than from the rest of the sky, this mod- 
ulating correction factor of diffuse radiation includes influence of 
cloudiness and it is expressed by: 


f= 2 (23) 


Ai, is the anisotropy index used to estimate the amount of cir- 
cumsolar diffuse radiation, also called forward scattered radiation. 
This factor is expressed as a ratio of beam radiation on a horizontal 
ground surface to extraterrestrial radiation. Its high value enhances 
the contribution of circumsolar diffuse radiation. The anisotropy 
index is given by: 


A= (24) 


2.4.3. Ground reflected radiation 
All models assume that the ground reflected component is 
isotropic as follows [49,57,58,75]: 


G=Gxpx (19) (25) 
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where p is the albedo coefficient of the ground. 


2.4.4, Global radiation 

The global radiation on a tilted surface, Gr, is the sum of the dif- 
fuse radiation, Gg, beam radiation Gp, and the ground reflected 
radiation, G,. Therefore, the incident global radiation on tilted sur- 
face is given by the following expression [76-81]: 


Gt = Got + Gar + Gr (26) 


In this study, the HDKR model will be applied assuming that the 
global radiation on tilted surface is: 


Gt = (Gp + G4A;)Rp + Ga(1 — A;)Rq rai’ (5) +GpgR- (27) 


3. PV generation 


The PV generation is mainly dependent on the global incident 
radiation (which was overviewed in the previous section) and the 
PV cell temperature. 


3.1. PV cell temperature 


It is well known that the PV cell temperature has an important 
effect on the PV output power. This temperature can be the same 
as the ambient temperature during the night but it can go over the 
ambient temperature by 30° or more in full sun. Hence, it is neces- 
sary to calculate the PV cell temperature in order to account for this 
meaning effect. For this reason, we have to begin by establishing the 
energy balance between, on one hand, the solar energy absorbed by 
the PV array, and on the other hand, the electrical output plus the 
heat transfer to the surroundings. The balance on unit area of the 
module which cooled by losses to the surroundings can be written 
as [33]: 


taGr = NcGr + UL(Tc — Ta) (28) 


where T, œ, Gr, nc, U, and Ta are respectively the solar transmittance 
of the PV array, the solar absorptance of the PV array, the global 
radiation striking the PV array, the electrical conversion efficiency 
of the PV array, the coefficient of heat transfer to the surroundings 
and the ambient temperature. 

According to the above equation, the PV cell temperature can be 
expressed as follows: 


ee Gr (3) (1 1c) (29) 


UL TO 


To estimate the value of (ta/U,), we report the nominal 
operating cell temperature (NOCT), which is defined as the cell 
temperature that results at an incident radiation of 0.8kW/m2, 
an ambient temperature of 20°C, and no load operation (meaning 
Nc =0), average wind speed equal to 1 m/s [82,83]. We can substi- 
tute these values into the above equation and solve it for (ta/U,): 
ta — Tc,noct — Ta,Nocr 


= 30 
UL GT, NOCT pau 


where Tcenocrt, Tanocr and Grnocr are respectively the nominal 
operating cell temperature, the ambient temperature and the solar 
radiation at which the NOCT is defined. 

If we presume that (tæ/U,) is constant, we can substitute this 
equation into the cell temperature equation as follows: 


T, -T, 
Te = Ta + Gr (‘ C,NOCT aoc) (1 w) (31) 


GT,NOCT TA 


We suppose that the PV array always operates at its maximum 
power point, which means that the cell efficiency is always equal 
to the maximum power point efficiency nmp: 


Nc = "mp (32) 


In this order, we can substitute nc by "mp in the relation (31) 
which can be written therefore as follows [84]: 


G 
Tc = Ta + (Tc,nocr — Ta,Nocr) (z z) (1 “np \ (33) 


However, nmp depends on the cell temperature Tc, we can 
assume that the efficiency varies linearly with temperature accord- 
ing to the following equation: 


Nmp = Nmp,stcl1 + &p(Tc — Tc,ste)] (34) 


where mp, œp and Tc stc are respectively the maximum power point 
efficiency under standard test conditions (a radiation of 1 kW/m2,a 
cell temperature of 25 °C, and no wind), the temperature coefficient 
of power and the cell temperature under standard test conditions. 
We note here that the temperature coefficient of power (œp) is neg- 
ative which means that the efficiency of the PV array decreases with 
increasing cell temperature. 

Finally, by substituting this efficiency equation (34) into 
the above cell temperature equation (33), the cell temperature 
becomes: 

Ta + (Tc,nocr — Ta,noct (Gs/Gs,nocr)[1 — (Nmp,stc/TÆ)(1 —JTc,stc)] (35) 


Ti 
i 1 + (Tc,nocr — Ta,nocr (Gs /Gs,noct )UMmp,stc / TŒ) 


3.2. PV output power 


The assessment of the PV output power was the interest of many 
authors. Thus, a several models [85] have been developed in order 
to evaluate all PV characteristics. In fact, Osterwald [86] suggested 
a method to calculate the power provided by a photovoltaic gener- 
ator. This method assumes that: 


Pm = Pm,stc (= ) [1 — æp(Tc — Teste) (36) 
Gr. stc 
where Pm and Pmstc are respectively the cell maximum power and 
the cell maximum power in the standard test conditions. 
Green and Araujo et al. [87,88] described another method to 
calculate the output power of a PV module as follows: 


RV kext a?) [1 2 x In(a)—rsx(1—a-)] (37) 


oc 
where Voc: cell open circuit voltage; Isc: cell short circuit cur- 
rent; oc: standard cell voltage; rs: standard resistance; a: 
Poc +1—2(Voc x rs); bi a/(1 +a). 

Araujo et al. [88] used the diode model to calculate the operation 
values of a PV generator. This model defines the V-I characteristic 
of the PV generator as follow: 


h exp [e — Voc + lg: “MeN | 


Ic = NmpNcplsc V: 


(38) 


where Ic: PV generator current (A); Vg: PV generator voltage (V); 
Ncp: parallel cells number of PV module; Ncs: parallel modules num- 
ber of PV generator; Nmp: parallel modules number of PV generator; 
Nms: series modules number of PV generator; R;: series resistance 
of a photovoltaic cell module (Q2); Vt: thermal voltage (V). 

In this study, we will use one of the simplest model to predict 
the PV array output power which assumes that: 


G 
Ppv = Ypvfev ( me 


) [1 + ep(Tc — Testc)] (39) 
GT, stc 
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Fig. 3. Instantaneous PV output power (f= 36°, y =0°) during the different equinox 
and solstice days. 


where Ypy, fev, Grstc, Œp, Tc and Tcstc are respectively the rated 
capacity of the PV array, meaning its power output under standard 
test conditions, the PV derating factor (it is used to account for such 
factors as shading, snow cover, aging, and so forth), the incident 
radiation at standard test conditions, the temperature coefficient 
of power (which indicates how strongly the PV array power out- 
put depends on the cell temperature), the PV cell temperature in 
the current time step and finally the PV cell temperature under 
standard test conditions. 


4. Outcome results 


This study is related to the Monastir city which is situated on the 
central coast of Tunisia, in the Sahel area (20 km south of Sousse and 
162 km south of Tunis). The latitude and longitude of this town are 
respectively 36°N and 12°E and an altitude of 10 m above sea level. 
Despite of meteorological data lack, we have used the same level 
of clearness index for all studied cases (K=1) basing on the Erbs 
et al. [54] clearness index correlation. It should be noted also that 
all simulations in this work are computed with a time step equal to 


((1/10) h). 
4.1. Effect of date on the instantaneous PV output power 


In order to study the dependence of the photovoltaic production 
with the effect of date throughout the whole year, four extreme 
days of the year are considered, i.e. vernal and autumnal equinox, 
and summer and winter solstice days. The results are carried out 
for a PV module oriented straight to the south by 0° Azimuth angle 
as well as a tilt angle equal to the latitude of Monastir (36°). 

The simulated results in terms of instantaneous PV output 
power for the abovementioned dates are illustrated in Fig. 3. Due to 
the difference in the length of the day between summer and winter, 
one can observe that the PV panel is generating electricity slightly 
for less than 10h at the winter solstice day and to slightly for more 
than 12 h at the summer solstice day. During the day, the PV output 
power will begin to increase from sunrise to midday then decline 
until sunset. Itis worth noting that the PV output power maximums 
registered at midday differ from date to another. This can be related 


to how high in the sky the sun appears. The operating hours of PV 
system in winter is shorter than in spring (or autumn) because the 
short sunshine time in winter compared to the vernal or autumn 
periods. This phenomenon was observed also by other researchers 
such as Chang [75] who found that in the northern hemisphere, 
the sun’s elevation angle during winter (November, December, and 
January) is lower than during the other seasons. This is because the 
mean distance between the Earth and the Sun in winter is less than 
that in spring (or autumn) by about 3% due to the eccentricity of 
the Earth’s orbit, resulting in a stronger flux of radiation. 


4.2. Effect of tilt angle on the instantaneous PV output power 


One of the important parameters that affect the performance 
of a PV system is its tilt angle with the horizontal. This is due to 
the fact that the variation of tilt angle changes the amount of solar 
radiation reaching the PV panel. 

Fig. 4 shows the instantaneous PV output power for different 
inclinations in Monastir during the four mentioned days (the winter 
solstice, the summer solstice, the vernal equinox and the autumnal 
equinox day). 

It is observed that the tilt angle has a very significant effect on 
the instantaneous photovoltaic production during the summer sol- 
stice day relatively compared to the effect during the winter solstice 
day, near the results of Hussein et al. [38]. This is because the zenith 
angle of the sun varies from 90° at the sunrise and the sunset to 
slightly more than 0° at the midday in the summer season while 
it varies only from 90° to slightly less than 60° in the winter sea- 
son. Thus, for a tilt angle equal 36°, the variation of solar radiation 
incident angle on the PV array is more important in summer days 
than in winter days. That’s why in the summer solstice day, the 
maximum PV output power varies from less than 60 W for an incli- 
nation angle equal to 90° to more than 140W for an inclination 
angle equal to 15° while in the winter solstice day, it varies only 
from more than 40 W for an inclination angle equal to 0° to less than 
90 W for an inclination angle equal to 60°. In the two equinox days, 
the maximum PV output power varies from more than 90 W for an 
inclination angle equal to 90° to less than 130 W for an inclination 
angle equal to 45°. 


4.3. Monthly optimal tilt angles for a south facing panel 


In order to enhance the PV output power of a module with- 
out tracking system, Gunerhan and Hepbasli [41] recommends that 
solar collector systems have to be adjusted with the optimal angle 
every month. Table 1 recapitulates the monthly optimal tilt angle 
for a PV module mounted facing the south in the Monastir city. The- 
ses angles vary from 4° in June to 65° in December. The variation 
of the optimal monthly tilt angle during the year is shown in Fig. 5. 
It is observed that in the first six months, the optimal tilt angle has 
the trend to decrease from 63° in January to 4° in June. However, in 
the second half of the year the optimal tilt angle witness a progress 
till 65° at the end of the year. This monthly optimal tilt angle ten- 
dency is in good agreement with the results of many researchers 
[33,41,53]. One can deduce that this variation has the same trend 
as the variation of the zenith angle of the sun during the year. In 
Monastir city and in other mid-latitude countries north of the equa- 
tor, the sun’s daily trip is an arc across the southern sky. The sun’s 
greatest height above the horizon occurs at noon, and how high 
the sun then gets depends on the season of the year, it is the high- 
est in mid-summer (June), which correspond to the lowest zenith 
angle value, the lowest in mid-winter (December) that correspond 
to the highest zenith angle value. This phenomenon was concluded 
also by Kelly and Gibson [89] in term of altitude angle. Indeed, they 
found that for the winter solstice day, the maximum altitude angle 
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Fig. 4. Instantaneous PV output power in terms of different inclinations during winter solstice (a), summer solstice (b), autumnal equinox (c) and vernal equinox (d) days. 


is @=24.2° while for the summer solstice day, the maximum solar Fig. 6 shows the yearly Output Energy at different tilt angles. 
altitude angle is œ =71.1°. For a horizontal PV Surface, the computed output energy is equal 
to 234kWh/year while for a vertical PV surface mounted to the 
south face the computed output energy is less than 167 kWh/year. 


4.4. Yearly tilt angle for a south facing panel It is found that the maximum of the yearly output energy can be 

obtained with a tilt angle equal to 33° which is close to 0.9 times 

Unfortunately, the cost of a monthly tracked PV System pro- latitude of Monastir city. This result is in good agreement with the 

hibits the use of this technology for many applications. This simulation result of Chang [58] who predicted that the yearly opti- 

inconvenience can be solved by finding out, with precision, the mal tilt angle, for a south-facing fixed panel, for latitudes below 65° 

yearly optimal tilt angle for a fixed PV module. in the northern hemisphere is equal to 0.9 latitude. For this angle, 
Table 1 


Monthly optimal tilt angles. 


Month January February March April May June July August September October November December 


Optimal slope 63° 55° 42° 25% 117 4 a 19° 34° 51° 61° 65 
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Fig. 5. Monthly optimal tilt angles for a facing south panel. 


the output energy can reach more than 268 kWh/year. However, 
for fixed PV using this value of tilt angle the energy produced in 
summer is not the maximum since it is calculated for the whole 
year and not for only the summer season. 


4.5. Comparison between the output power generated by a 
traditional fixed and single-axis tracked panel 


A comparison of monthly output energy between fixed and 
tracked southern face mounted PV module is illustrated in Fig. 7. 
It is noticed that the optimal yearly tilt angle computed below 
(i.e. 33°) is used for the fixed panel and the optimal monthly 
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Fig. 6. Yearly optimal tilt angle for a facing south panel. 
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Fig. 7. Monthly PV output power generated by a traditional fixed and a single-axis 
tracked panel. 


tilt angles are used for the tracked system. One can observe 
that the tracked panel generates stronger output energy than a 
traditionally fixed panel due to its smaller incidence angle of sun- 
light. 

The monthly increment of output energy obtained by the 
tracked panel is marked with black color. One can observe that the 
gain provided by the tracked panel is more noticeable for tow crit- 
ical periods. The first period corresponds to the months of May, 
June and July and the second period corresponds to November, 
December and January months. The first period correspond to the 
summer solstice at June, in which the gain reveals its highest value 
equal to 2.9kWh corresponding to 10.34% of energy generated 
by a fixed PV collector. The second period marked by the occur- 
rence of the winter solstice at December. During this period, the 
gains made by the tracked panel relative to a fixed panel reach, 
in December, more than 15% corresponding to 2.06 kWh. Through- 
out the other months, the gain reaches only a value slightly more 
than 1kWh in the month of February. Table 2 lists the Output 
power of PV and gains for different time periods. This variation 
of Gain can be explained by the fact that in summer and winter 
solstices the solar zenith angles reach respectively their lowest 
and highest values so the incidence angle of solar rays on fixed 
panel is more important at these periods. The yearly accumulation 
gain made by the tracked panel relative to a fixed panel is eval- 
uated to 15.48kWh which correspond to 5.76%. This result was 
concluded by many other researchers. For instance, Lubitz [49] 
found that the gain made by tilt PV tracker across the contigu- 
ous United States of America with a value of 5% Neville. Tomson 
[90] described mainly the performance of PV modules with daily 
two-position in the morning and in the afternoon. His results indi- 
cated that the seasonal energy yield was increased by 10-20% over 
the yield from a fixed south facing collector tilted at an optimal 
angle. 


4.6. Effect of azimuth angle on the instantaneous output power 


When a module panel must be installed with the azimuth other 
than south, Fig. 8 illustrates the instantaneous output power for 
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Table 2 
Output power of PV and gains for different time periods (kWh). 


Month January February March April May June July August September October November December Total 

Fixed panel (kWh) 14.64 16.91 23.76 26.6 28.85 28.04 28.85 27.98 24.31 20.27 15.05 13.34 268.60 
Tracked panel (kWh) 16.52 17.99 23.97 26.8 30.55 30.94 31.16 28.60 24.32 21.10 16.73 15.40 284.08 
Gain (kWh) 1.88 1.08 0.21 0.20 1.70 2.90 2.31 062 0.01 0.83 1.68 2.06 15.48 
Gain (%) 12.84 6.38 0.88 0.75 5.89 10.34 8.00 2.21 0.04 4.09 11.16 15.44 5.76 


different azimuths respectively during summer and winter solstice 
days and autumnal and vernal equinox days. The panel, here, is 
installed with a constant tilt angle equal to the yearly optimal tilt 
angle (i.e. 33°) for the city of Monastir. In order to analyze the 
azimuth effect, many azimuth angles are investigated, i.e. —120°, 
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—80°, —40°, 0°, 40°, 80°, and 120°. We remember that the azimuth 
angle due to the south is 0°, due east is —90°, due west is 90°, 
and due north is 180°. With fixed-azimuth systems, the panels 
are almost oriented towards the equator (0° azimuth in the north- 
ern hemisphere). One can observe that varying the azimuth angle 
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Fig. 8. Instantaneous PV output power in terms of different azimuth angles during winter solstice (a), summer solstice (b), autumnal equinox (c) and vernal equinox (d) days. 
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from —120° to 120° modify the position and the value reached by 
the maximum output power. For azimuth angle equal to 0° the 
peak of generated power is registered at the solar noon and it is 
the highest value relatively compared to the values reached for 
other azimuths. For negative values of the azimuth angles (east) 
the output power maximums are reached before the solar noon 
while for positive values of the azimuth angles (west) the output 
power maximums are reached after the solar noon. In summer 
solstice day, the maximum output power is recorded at 10h54’, 
11h00’, 11h30’, 12h12’, 13h00’, 13h24’ and 13h30’ respectively for 
azimuth angles equal to —120°, —80°, —40°, 0°, 40°, 80° and 120°. 
Theses translations in time are accompanied by a decrease in term 
of the maximum output power values. In summer solstice day, this 
decrease is significant only for azimuth angles equal to —120° and 
120° while in winter solstice day the decrease is significant for all 
considered azimuth angles. The decrease do not exceed 8.03% for 
azimuth angle equal to —120° or 120°, it is evaluated to 1.45% for 
azimuth angle equal to —80° or 80° and it is less than 0.01% for 
azimuth angle equal to —40° or 40° in summer solstice day while 
in winter solstice day this reduction achieve more than 65%, 32%, 
7% respectively for azimuth angle equal to +120°, to +80°, to +40°. 
For vernal and autumnal equinox days, the decrease is between 
3.9% and 35% for azimuth angle between +120° and +40°. These 
results do not ignore the importance of installing PV panel facing 
a direction other than the south specially to cover a load at earlier 
or later hour. One can observe that for example in summer solstice 
day, the instantaneous output power for a PV panel with azimuth 
angle equal to 0° is evaluated only to 14W and 4 W respectively at 
7h and 18h and it can reach the values of 52 W and 34W respec- 
tively at 7h and 18h with azimuth angles equal respectively to 
—120° and 120°. Hence, the gain made by choosing theses azimuth 
angles compared to the south direction is estimated to 271% and 
750%. 


4.7. Investigation of optimal tilt and azimuth angles 


In the above analysis, we investigated separately the effect of tilt 
and azimuth angles on the instantaneous PV panel output power 
in the city of Monastir. In order to optimize the exploitation of PV 
installations, dual-axis tracked systems can be used. The first axe is 
a horizontal one to adjust the tilt angle, the second axe is a vertical 
one to adjust the azimuth angle. So the determination of a couple 
of optimal angles during the day has to be examined. Fig. 9 illus- 
trates the variation of optimal tilt and azimuth angles of PV panel, 
in Monastir city, during four extreme days: winter and summer 
solstice days and autumnal and vernal equinox days. The results 
show that the optimal tilt angle is equal to 90° at the sunrise and 
at the sunset. A tilt angle equal to 90° corresponds to a vertical PV 
panel oriented to the sun at the horizon. In the first half of the day, 
this angle decrease from 90° until reaching the minimum at the 
solar noon, in the second half of the day it increases until 90° at 
the end of the day. This phenomenon was observed also by Chang 
[58] who showed that the optimal tilt angles of the tracked panel 
is O° at noon, and increases towards dawn and dusk. The minimum 
reached by the optimal tilt angle at the solar noon varies from a 
day to another. Its highest value (i.e. 59°) is obtained for the winter 
solstice day and its lowest value (12°) is obtained for the summer 
solstice day. At any hour during the day, the optimal tilt angle for 
the summer solstice day is lower than the one for the winter sol- 
stice day. This is because in the winter season the sun zenith angle is 
more important than in the summer season. So the PV panel should 
be more inclined to be perpendicular to the sun rays. Besides, a tilt 
angle equal to 90° is suggested at 8 h and 17 hin the winter solstice 
day while this inclination is recommended before 7 h and after 19h 
in the summer solstice day. This is related to the above stated fact 
that the day is longer in summer than in winter season. 
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Fig. 9. Instantaneous optimal tilt and azimuth angles during winter solstice, sum- 
mer solstice, autumnal equinox and vernal equinox days. 


During the day, the optimal azimuth angle increases from neg- 
atives values corresponding to the east direction at the morning to 
positives values corresponding to the west direction. In the winter 
solstice day, it varies only from —70° to 70° wile in the summer sol- 
stice day it varies from less than —110° to more than 110°. This is 
explained by the same reasons mentioned above. One can observe 
that the variation rate of the optimal azimuth angles is more impor- 
tant between 10h and 14h. At this period surrounding the solar 
noon, the solar trajectory in the sky is perpendicular to the axe 
joining the PV panel to the sun. Thus the regular rotation of the 
earth (15°/h) corresponds to the highest displacement of the sun 
on its trajectory in the sky. This large displacement recommends a 
large variation of the azimuth angle. 


4.8. Comparison between the output power generated by a 
traditional fixed and dual-axis tracked panel 


Fig. 10 shows the instantaneous PV output power generated by 
fixed module with the yearly optimal tilt angle (i.e. 33°) facing due 
to the south and a dual axis tracked panel during the winter solstice 
(a), summer solstice (b), autumnal equinox (c) and vernal equinox 
(d) days, respectively, in Monastir city. It is found that the two- 
axis tracked panel produces stronger power than a traditionally 
fixed panel due to its optimal incidence angle of sunlight during 
the day. The same conclusion was proposed by Kacira et al. [91] 
who investigated the performance of PV modules mounted on a 
two-axis solar tracker. They concluded that the gain in the daily 
power generation could be 34.6%. The instantaneous increment of 
output power obtained by the dual-axis tracked panel is marked 
with black color. A similar gain was concluded by Tina and Gagliano 
using a two-axis tracked system in the site of Catania (Italy) [92]. 
The increment in term of PV produced energy made by the dual 
axis tracked PV system during the winter and summer solstice days 
reaches respectively 128 Wh and 413 Wh. Theses increments made 
by the dual-axis tracked PV system correspond to a gain of 30% and 
44% respectively in the winter and summer solstice days relatively 
compared to a traditionally fixed panel. 
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Fig. 10. Instantaneous PV output power generated by a traditional fixed and a dual-axis tracked panel during winter solstice (a) summer solstice (b), autumnal equinox (c) 


and vernal equinox (d) days. 


5. Conclusion 


An overview on research works on solar radiation basics 
and photovoltaic generation has been presented. Indeed, vari- 
ous models and approaches have been discussed to assess the 
total solar radiation on horizontal and tilt surface. Several meth- 
ods for evaluating the output power of a PV module have been 
presented. A complete performance modeling and investigation 
of fixed, single and dual-axis tracking photovoltaic panel have 
been evaluated with details. These formulations are very useful 
for PV system assessment and planning for PV projects. Using 
the case study of Monastir, some concluding remarks are as 
follows: 


e For Monastir city, the yearly optimal tilt angle of a fixed panel 
faced due to the south is close to 0.9 times the latitude. 


e The gain made by the module mounted on a single-axis tracking 
panel relative to a traditional fixed panel are more noticeable 
for two critical periods which correspond to those surrounding 
the summer and the winter solstice dates. It reaches the value 
of 10.34% and 15% in the summer and winter solstice periods, 
respectively. 
The yearly gain made by a module mounted on a single-axis 
tracking system relative to a fixed panel installed with the yearly 
optimal tilt angle is estimated to 5.76%. 
In order to more optimize the solar systems exploitation, a dual- 
axis tracking installation can be used. The gains made by a 
dual-axis tracking panel relative to a traditional fixed panel are 
investigated. This gain reaches 30% and 44% respectively in the 
winter and summer solstice days. 
e This study provide a guide to the choice of the optimal tilt and 
azimuths angles for fixed and tracked panel mounted on sin- 
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gle and dual-axis photovoltaic systems used in Monastir city 
(Tunisia). 
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